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ABSTRACT 



We numerically investigate orbital evolution of star clusters (SCs) under the influ- 
ence of dynamical friction by field stars of their host disk galaxies embedded in dark 
matter halos. We find that SCs with masses larger than ~2x 1O 5 M can show signif- 
icant orbital decay within less than 1 Gyr due to dynamical friction by disk field stars 
in galaxies with disk masses (Mi) less than 1O 9 M . We also find that orbital decay of 
SCs due to dynamical friction is more remarkable in disk galaxies with smaller Md and 
higher mass-ratios of disks to dark matter halos. The half-number radii (i?h,sc) and 
mean masses within i?h,sc of the SC systems (SCSs) in low-mass disk galaxies with 
Md ^ 1O 9 M0 are found to evolve significantly with time owing to dynamical friction 
of SCs. More massive SCs that can spiral-in to the central regions of disks can form 
multiple SC systems with smaller velocity dispersions so that they can merge with 
one another to form single stellar nuclei with their masses comparable to ~ 0.4% of 
their host disk masses. Based on these results, we suggest that luminosity functions 
(LFs) for more massive globular clusters (GCs) with masses larger than 2 x 10 5 Mq 
can steepen owing to transformation of the more massive GCs into single stellar nuclei 
through GC merging in less luminous galaxies. We also suggest that the half-number 
radii of GC systems can evolve owing to dynamical friction only for galaxies with their 
total masses smaller than ~ 1O 1O M0. 

Key words: globular clusters:general - galaxies:formation - galaxies:kinematics and 
dynamics - galaxies: halos - galaxies:star clusters 



1 INTRODUCTION 

Dynamical friction of SCs in galaxies have been discussed in 
many different contexts, such as formation of stellar galac- 
tic nuclei via merging of old GCs (e.g., Tremaine et al. 
1975), transformation from non- nucleated dwarfs into nu- 
cleated ones (e.g., Oh & Lin 2000), physical meanings for 
the presence of the GC system (GCS) in the Fornax dwarf 
galaxy (e.g., Oh et al. 2001), and dark matter distributions 
of dwarf galaxies (e.g., Hernandez & Gilmore 1998; Goerdt 
et al. 2006; Inoue 2009). Orbital decay of GCs due to dy- 
namical friction can significantly change spatial distribution 
of GCs and thus increase possibilities of them to be de- 
stroyed by strong galactic tidal fields (e.g., Vesperini 2000). 
Dynamical friction of SCs (including GCs) in galaxies is 
thus considered to be important for better understanding 
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the evolution of GC luminosity functions (GCLFs) and that 
of half- number radii of the GCSs in galaxies (e.g., Vesperini 
2000). Recently, van de Ven & Chang (2009) have suggested 
the importance of dynamical friction in the dynamical evo- 
lution of SCs around nuclear rings in galaxies. 

Dynamical friction processes of SCs against disk field 
stars in disk galaxies are important for the following three 
reasons. Firstly, recent cosmological hydrodynamical simu- 
lations on the possible formation sites of first GCs (Kravtsov 
& Gnedin 2005) have shown that the present-day GCs can 
be formed within disk galaxies at high redshifts (2 ~ 3). Sec- 
ondly, less luminous disk galaxies like the Large Magellanic 
Cloud (LMC) have disky GC systems in which GCs have 
disky spatial distributions and rotational kinematics (e.g., 
Freeman et al. 1983; Olsen et al. 2004). Thirdly, previous 
theoretical studies suggested that the present-day dwarf el- 
lipticals (dEs) were previously less luminous disk galaxies 
with no/little bulges (Mastropietro et al. 2004). Further- 
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Table 1. The ranges of model parameters. 

Parameters M ga i (xlO lo M )» M d (xlO 9 A/ )^ f d ^ f b £ Galaxy type ^ Af v ,iowX M VMe h£. 

Value ranges 0.1-10.0 0.1-10.0 0.01-"o.l 0.0-L0 LSB or HSB™ -6 (mag)" -12 ~ -8 (mag) 

a The total mass of a galaxy. 
^ The total mass of a disk. 

c The mass fraction of stellar disk (= M^/M sa \) in a galaxy. 
^ The mass fraction of a bulge (= M^/Md) in a galaxy . 

e HSB and LSB represent low-surface brightness and high-surface brightness stellar disks, respectively. 

J The lower luminosity cut-off in the SC luminosity function. This parameter is fixed at —6 mag in all models (see the main text for 
reasons for this). 

9 The higher luminosity cut-off in the SC luminosity function. 



more, previous observational results on physical properties 
of dwarf galaxies (e.g., Stiavelli et al. 2001) suggest simi- 
larity in nuclear density profiles of stars between dEs and 
spirals with exponential bulges: these appear to imply that 
understanding nucleus formation processes due to SC mi- 
gration caused by dynamical friction in spirals can further 
lead to better understanding those in dwarfs. Although dy- 
namical friction processes of SCs against dark matter halos 
have been well investigated by previous works (e.g., Goerdt 
et al. 2006), those by field stars in disk galaxies have not 
been done by self-consistent numerical simulations so far. 

The purpose of this paper is thus to investigate dynam- 
ical friction of SCs against disk field stars based on fully 
self-consistent numerical simulations on dynamical evolution 
of disk galaxies with SCs. We particularly investigate how 
dynamical friction of SCs against disk field stars depend 
on physical properties of their host galaxies, such as disk 
masses, bulge-to-disk-ratios, and mass-ratios of disk to dark 
matter halo. Based on the present numerical results, we dis- 
cuss the observed luminosity functions of GCs dependent on 
luminosities of their host galaxies (e.g., Jordan et al. 2006), 
formation processes of stellar galactic nuclei, and evolution 
of physical properties of GCSs in galaxies. 

The plan of the paper is as follows: In the next section, 
we describe our numerical models for orbital evolution of 
SCs in disk galaxies. In §3, we present the numerical results 
mainly on (i) orbital evolution of SCs and (ii) the physi- 
cal properties of SCSs for variously different models. In §4, 
we discuss wide implications of the present results such as 
formation of stellar galactic nuclei by SC merging and evolu- 
tion of mass and luminosity functions of SCs due to orbital 
decay of SCs caused by dynamical friction. We summarize 
our conclusions in §5. 



2 THE MODEL 

We investigate orbital evolution of SCs in disk galaxies em- 
bedded in massive dark matter halos using the latest version 
of GRAPE (GRavity PipE, GRAPE-7) which is the special- 
purpose computer for gravitational dynamics (Sugimoto et 
al. 1990). All SCs are represented by point-mass particles 
rather than self-gravitating N-body systems so that we can 
not discuss internal stellar dynamics of SCs influenced by 
galactic potentials. Although evolution of SCs (e.g., internal 
stellar dynamics) from their formation within GMCs in a 
low-mass disk galaxy is discussed by our previous numerical 
simulations (e.g., Hurley & Bekki 2008), the present numer- 



ical code does not allow us to discuss this important issue: 
we plan to discuss this in our forthcoming papers. 

Galaxies are represented by N-body particles so that 
SCs can feel gravitational influences of the live galactic po- 
tentials (e.g., dynamical friction). Non-axisymmetric struc- 
tures such as bars and spiral arms in disk galaxies can change 
background stellar distribution and kinematics, which can 
be important for the effectiveness of dynamical friction 
that depends on density and velocities of background stars. 
Therefore, dynamical friction of SCs against disk field stars 
can be very complicated in the present study. Dynamical 
friction of disk SCs (i.e., SCs initially in stellar disks) against 
disk field stars is much more effective than that against dark 
matter halos in the present study. Therefore, the orbital de- 
cay of disk SCs demonstrated in the present study is due 
mainly to gravitational interaction between disk field stars 
and SCs. 

2.1 Disk galaxy model 

Since our numerical methods for modeling dynamical evolu- 
tion of late- type disk galaxies have already been described by 
Bekki & Peng (2006), we give only a brief review here. The 
total disk mass and the size of a disk of a disk galaxy with 
the total mass of M ga i are Md and Rd, respectively. Time 
is measured in units of t g = (i?j/GM ga i) 1//2 , where G is the 
gravitational constant and assumed to be 1.0. If we adopt 
Aigai = 1.0 x 10 10 M Q and R d = 2.3 kpc as a fiducial value, 
then t g = 1.6 x 10 7 yr. The model with this fiducial value is 
referred to as "the standard model". The disk is composed 
of a dark matter halo, a stellar disk, and a stellar bulge. We 
mainly investigate late-type, bulge-less galaxies, because we 
mainly discuss SC evolution in low-mass galaxies. 

The mass ratio of the dark matter halo (with the mass 
of Md m ) to the stellar disk in a disk model is a free parameter 
ranging from 9 to 49.0 in order to investigate galaxies with 
different M ga i thus different mass-to-light-ratio. We adopt 
the density distribution of the NFW halo (Navarro, Frenk 
& White 1996) suggested from CDM simulations: 

PW = (r/r s )(l + r/r s )2' (1) 

where r, po, an d r B axe the spherical radius, the characteristic 
density of a dark halo, and the scale length of the halo, 
respectively. The c parameter (= r s /r v i r , where r v i r is the 
virial radius of the NFW profile) for a galaxy with Mdm 
is chosen according to the predicted c-Md m relation in the 
ACDM simulations (e.g., Neto et al. 2007). For example, 
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reasonable c values for disk galaxies with M ga i = 10 10 Mq 
are 12.9 (also r v j r is ~ 10J?d)- The stellar bulge with a mass 
Mb and size Rb is represented by the Hernquist profile with 
the scale-length of 0.2i?b- The bulge mass fraction (Mb/Ma) 
is referred to as /b- 

The radial (R) and vertical (Z) density profiles of the 
disk are assumed to be proportional to exp(— R/Ro) with 
scale length Ro = 0.2 and to sech 2 (Z/Zo) with scale length 
Zq = 0.04 in our units, respectively: the stellar disk follows 
this exponential distribution. In addition to the rotational 
velocity caused by the gravitational field of disk, bulge, and 
dark halo components, the initial radial and azimuthal veloc- 
ity dispersions are assigned to the disc component according 
to the epicyclic theory with Toomre's parameter Q — 1.5. 
The vertical velocity dispersion at given radius is set to be 
0.5 times as large as the radial velocity dispersion at that 
point, as is consistent with the observed trend of the Milky 
Way (e.g., Wielen 1977). 

We investigate models with different M ga j and adopt 
the Freeman's law (Freeman 1970) to determine Ro of a 
disk galaxy according to its disk mass: 



R 



■■ 3.5(- 



M d 



0.5 

-) kpc. 



(2) 



6 x 1O 1O M ' 

Structural and kinematical properties of dark matter ha- 
los and stellar disks are assumed to be self-similar between 
models with different M ga i. We also investigate low-surface 
brightness (LSB) models with different M ga i in which Ro of 
a disk galaxy is determined as: 

M d 



Ro = 8.. 



0.5 

-) kpc. 



(3) 



v 6 x 10 10 M Q 

Both the total particle number of dark matter halo (iVh) 
and that of the stellar disk (N c \) in a bulgeless disk model are 
500000. The total number for bulge (iVb) in a model depends 
on Mb such that Nb = iVa Mb/Ma- The fixed gravitational 
softening length is 0.008i?<j, which corresponds to 19pc in 
the standard model with Md = 1O 9 M0. The mass-ratio of 
the lowest-mass SC to each disk field star is more than 10 
owing to the adopted particle number (~ 10 6 ) so that we 
can properly investigate dynamical friction processes of SCs 
by disk field stars. The leapfrog integration scheme with a 
fixed time step of 0.02f g is adopted for all models. 

2.2 SC model 

Since we focus on orbital evolution of long-lived SCs in disk 
galaxies, we adopt the canonical luminosity function (LF) 
observed for old GCs (Harris 1991): 



$(M) = C x 



- (M - M ) 2 /2<j„ 



(4) 



where C is a constant, M is the magnitude of a SC, 
M =-7.27 mag in the V band (Harris 1991), and cr m = 1.2 
mag. We assume that Mo and <r m should be the same be- 
tween different models in order to show more clearly how 
dynamical friction processes of SCs depend on physical prop- 
erties of their host galaxies. The mass function (MF) of SCs 
(SCMFs) is derived from the LF os SCs (SCLFs) for a rea- 
sonable stellar- mass-to-light-ratio. 

The V-band magnitudes of SCs for lower and higher 
cut-off luminosities/masses in SCLFs/SCMFs are described 
as Mv.iow and My, high, respectively. Mv,iow and Mv.high are 
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Figure 1. Spatial distributions of disk SCs represented by small 
black dots in the disk galaxy projected onto the x-y plane for 
four different time steps in the standard model. Upper and lower 
panels show the entire and nuclear distributions of SCs, respec- 
tively. Time T in units of Gyr is given in each of the four panels. 
Big blue and small red circles represent the initial disk size of the 
galaxy and the nuclear region (Ft = 100 pc), respectively. 



set to be —6 mag and —10 mag, respectively for most mod- 
els. Although we consider these values to be reasonable and 
realistic, we investigate models with different My, high for 
comparison. If we adopt Myjow significantly larger than —6 
mag, the mass of a SC can be comparable to each disk field 
star in some models so that we can not properly investigate 
dynamical friction processes. We therefore consider that it is 
the best for the present study to fix Mv,iow at a reasonable 
value. The stellar-to-mass-to-light-ratio is 3.2 for SCs in all 
models to allocate a mass (m sc ) to each SC particle. 

The total number of SCs in a galaxy is determined by 
Sn of the galaxy, where SW is the specific frequency of the SC 
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Figure 2. Orbital evolution of a SC with m sc = 3.3 X 10 5 M Q 
projected onto the x-y plane (left) and the x-z plane (right). The 
initial and final locations of the SC are indicated by a blue square 
and by a red circle, respectively. The orbit for ^ T ^ 20t g is 
shown so that the orbital decay due to dynamical friction can be 
more clearly seen. 
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Figure 3. The same as Fig. 2 but for m BC = 1.3 X 1O 6 M . 



smaller/larger than those consistent with the adopted 5n- 
My relation for comparison. 

We assume that the initial distribution and kinematics 
of the SCS in a disk galaxy are exactly the same as those of 
the disk components (i.e., rotating, exponential disk). This 
is a reasonable assumption, given that the vast majority of 
stars are formed as bound and unbound SCs in the Galac- 
tic disk (e.g., Lada & Lada 2003). For the adopted expo- 
nential distribution of a SCS, we allocate each SC particle 
3D-velocities as follows. We first search for the nearest disk 
particle (i.e., disk field star) for each SC particle and then 
allocate the 3D velocities of the disk particle to the SC one. 
Thus the initial SCS in a galaxy has a thin disk configuration 
and rotational kinematics. 

In order to discuss the effectiveness of dynamical fric- 
tion of disk SCs against disk field stars, we also briefly in- 
vestigate orbital evolution of halo SCs that are initially in 
halo regions of galaxies. The number of halo SCs is the same 
as that of disk SCs in a galaxy and its spatial distribution is 
described as a power-law with the power-law index of —3.5 
as observed for the Galactic GCs (e.g., van den Bergh 2000). 
The initial half-number radius of the halo SCs in a galaxy 
is 1.4Ro as observed for the Galactic GCS (i.e., the half- 
number radius of 5 kpc for the GCS and Ro = 3.5 kpc for 
the stellar disk; van der Bergh 2000). Kinematical proper- 
ties of the halo SCs are assumed to be the same as those of 
the dark matter halos of their host galaxies (i.e., isotropic 
velocity dispersions determined by mass distributions of the 
galaxies) . 



system. The specific frequency is defined as follows (Harris 
& van den Bergh 1981): 

5 N = iV sc x 10°- 4(M ^ +15) , (5) 

where N sc and M v are the total number of SCs in a galaxy 
and V— band absolute magnitude of the galaxy, respectively. 
Thus, Nac is determined by Sn and Md with a reasonable as- 
sumption on the mass-to-light-ratio for stars (for converting 
Md into M v in galaxies). 

The observed dependence of Sn on My has "U-shape" 
in the sence that the galaxy luminosity dependence of Sn 
is different between galaxies below and above a threshold 
luminosity, Afv,th, which is around —19.5 mag (Bekki et al. 
2006). We here use the following form for the Ss-My relation 
(Bekki et al 2006); 

Ss{x) =AxX W Klx + (S N , th - Ax) x 10 A ' 2 ", (6) 
where x — M \, Mv - th and Sn t h is .Sim at My = My th . 

M V,th 

Parameter values of A\, Ki, K2, Afy.th, and 5*N,th can be 
determined by fitting to observations. A reasonable model fit 
to observational data sets have Mv,th = —19.5 mag, SW.th = 
1.0, Ax = 0.5, Ki = -6, and K% = 4 (see Fig. 1 in Bekki et 
al. 2006). 

We mainly investigate the "standard" SC model in 
which M d = 10 9 M Q , M A /L = 4.0 (where L is the to- 
tal luminosity of the disk), Sn = 6.8 (thus N BC = 20). 
This model can be consistent with the observed U-shaped 
SW-Mv relation. The models with the same Ss-My rela- 
tion as that adopted in the standard model have iV sc = 3 
for M d = 1O 8 M and N sc = 117 for M d = 10 10 M Q in 
the present study. We however investigate models with N sc 



2.3 Parameter study 

We mainly investigate dynamical friction processes of SCs 
with different initial masses, their orbital evolution and the 
final distribution of SCs in the standard model with Md = 
1O 9 M , / b = 0, / d = 0.1, iV sc = 20, and Af v , h i gh = -10 
mag. We also investigate 40 models with different Md , /b , fd , 
N sc , Afy.high, central densities of stellar disk (i.e., whether 
disks are LSB or HSB, where LSB and HSB represent low- 
surface brightness and high-surface brightness disk galaxies, 
respectively) and structures of dark matter halos in order to 
better understand the effectiveness of dynamical friction in 
disk galaxies with different physical properties. The range 
of model parameters investigated are shown in the Table 1. 

In order to discuss formation of stellar nuclei in disk 
galaxies, we investigate time evolution of the total mass of 
SCs within O.liid for each model. We consider that SCs 
transfered to the central 0.1-Rd owing to dynamical friction 
can quickly merge with one another to form a single massive 
SC that can be identified either as a nuclear SC or a stellar 
galactic nucleus. Since SCs are represented by point-mass 
particles in the present simulations, nucleus formation via 
SC merging can not be investigated in a fully self-consistent 
manner. Previous and recent numerical simulations in which 
SCs are represented by N-body particles have confirmed that 
nuclear SCs can merge to form a single nucleus in a galaxy 
(e.g., Capuzzo-Dolcetta & Miocchi 2008). We thus consider 
that it is appropriate for the present study to briefly discuss 
formation of stellar galactic nuclei via SC merging in central 
regions of galaxies based on our simulations. 

Thus the total mass of the stellar nucleus in a disk 
galaxy at each time step in each simulation (M nuc ) is de- 
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scribed as follows: 



1.2 



M nuc = 



=1 



(7) 



where N nuc is the total number of nuclear SCs within 0.1i?d 
and m SCt i is the mass of the individual nuclear SCs. M nuc is 
not literally the nucleus mass of the galaxy so that it can 
change with time relatively rapidly. For all models, we run 
models for at least 40i g in order to investigate orbital decay 
of SCs due to dynamical friction in galaxies. We however 
investigate some models, in particular, more massive disk 
ones, for longer time scales (up to 200t g ) so that we can 
confirm whether dynamical friction is really ineffective in 
these massive disk models. 

The final M nuc in a galaxy depends strongly on the Sn 
such that Af nuc can be larger in models with larger Sn for 
a fixed Md- Although this suggests that the observe higher 
S N in nucleated dwarfs (dE,Ns; Miller & Lotz 2007) is due 
to initially high Sn in dE,Ns, we do not discuss extensively 
this problem based on results of models with different Sn 
(by adopting Sn different from those used in the observed 
Sn — My relation) : we will discuss this problem in our forth- 
coming papers. In the followings, the time T represents the 
time that has elapsed since the simulation starts. 



2.4 Comparison with analytical works 

Previous analytical works showed that dynamical friction 
processes of SCs within galaxies depend on masses of SCs 
and physical properties of their host galaxies (e.g. Binney 
& Tremaine 1987). These analytical works would be quite 
useful and helpful in better understanding numerical results, 
though they did not investigate dynamical friction of SCs 
against disk field stars. The time-scale of dynamical friction 
(tfric) for typical GCs with masses (M gc ) of 2xlO 5 M within 
galactic halos of luminous galaxies is considered to be quite 
long (e.g., Binney & Tremaine 1987). For the Galaxy with 
the circular velocity v c = 220 km s _1 and the Coulomb 
logarithm In A = 10, tf r i C can be estimated as follows: 



fa 



1.2 x lO 11 ^ 



;) (J 



-)(; 



") J*, (8) 



2kpc' '220kms- 1 "2 x 10 5 M Q ' 

where n is the distance of a GC from the center of its host 
galaxy. For a galaxy with v c — 70 km s _1 and In A = 10, 



iMc = 1.8 x 10 9 (- 



2 

:) (- 



-)(■ 



yr. 



(9) 



'lkpc y v 7Okms- iyv lQ 6 M ' 

Therefore, although dynamical friction can not be impor- 
tant for evolution of normal SCs and GCs within luminous 
galaxies like the Galaxy, it can significantly change physi- 
cal properties of SCs and those of SCSs in dwarf galaxies. 
The above results also show that dynamical friction of SCs 
can be more effective for SCs initially located in the inner 
regions of galaxies. 

These estimation is for dynamical friction of SCs against 
galactic halos so that we can not rely on these for under- 
standing the present numerical results for those against disk 
field stars. However, the above analytical estimation would 
be still useful for (i) understanding the differences in if r i c 
between dynamical friction against galactic halos and that 
against disk field stars and (ii) discussing dependences of the 
present results on Md (thus on « c ) in a qualitative manner. 
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Figure 4. The ratios of final R (R{) to initial one (Ri) as a 
function of ij; for SCs in the standard model. Blue and red dots 
represent SCs with m sc < 2 x 10 5 Mq and m sc > 2 X 10 5 Mq, 
respectively. The horizontal dotted line indicates the mean value 
of Rf/Ri for all SCs in this model. If we estimate the mean for 
SCs with Rf/Ri 1 for which orbital decay due to dynamical 
friction is clearly seen, the mean value is 0.55. 
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Figure 5. The ratios of minimum R (R mln ) to initial one (Ri) as 
a function of Ri for SCs in the standard model. Blue and red dots 
represent SCs with m sc < 2 x 1O 5 M0 and m sc ^ 2 X 1O 5 M0, re- 
spectively. The horizontal dotted line indicates the mean R m i n /Ri 
in this model. 



3 RESULTS 

3.1 The standard model 

Fig. 1 shows the time evolution of the projected distribu- 
tions of SCs for the last 0.64 Gyr in the standard model 
with Md = 1O 9 M and with iV sc = 20. The projected spatial 
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Table 2. A brief summary of the results for the five representative models with different disk masses. 



(xlO 9 M ),^ 


M nuc (xl0 6 M Q )^ 


M nuc /M d- ^ 




comments 


0.1 


0.5 


0.0049 


0.28 


low-mass disk model 


0.5 


2.2 


0.0043 


0.50 




1.0 


4.3 


0.0043 


0.55 


standard model 


5.0 


3.0 


0.0006 


0.67 




10.0 


2.8 


0.00028 


0.83 


high-mass disk model 



a The total mass of a disk. For these five disk models, other model parameters such as f d and /j, are exactly the same between the five. 
" The total mass of the stellar nucleus in a disk: This is the total mass of SCs within O.liZd- 
c The nuclear mass fraction of a disk. 

^ The mean value of Rj/R^, where i?; and Rf are initial and final radii of SCs from the center of their host galaxy. This mean value is estimated 
for SCs with Rf/R\ 1 so that we can more appropriately measure the degree of orbital decay of SCs due to dynamical friction. 
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Figure 6. Time evolution of the mass-ratio of the nucleus to the 
disk (M n uc/M(j) in the standard model. The dotted line is the 
observed mean of M nuc / 'M d , See the main text for the details of 
the method to estimate Afnuc 



0.8 - 



\ 0.6 - 

k : 



0.4 
0.2 



% 0.6 • • 

k : 



0.4 - 



0.2 



M d =5*10 B M s 



0.2 0.4 0.6 O.f 



1 0.2 0.4 0.6 0.8 1 

R, (*R„) 



Figure 7. The ratios of final R (Rf) to initial one (ii;) as a 
function of ij; for SCs in four models with M d = 10 s Mq (upper 
left), M d = 5 X 1O 8 M (upper right), M d = 5 X 1O 9 M (lower 
left), and M d = 1O 1O M (lower right). 
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Figure 8. The same as Fig. 7. but for four models with 
M dm /M d = 19 (upper left), M dm /M d = 49 (upper right), 
M b /M d = 0.17 (lower left), and M b /M d = 1 (lower right). 



distribution of SCs can change with time not only because 
of the dynamical evolution of the disk but also because of 
dynamical friction of SCs against the disk field stars. The 
half-number radius of the SCS evolves from 0.38J?d (cor- 
responding to = 0.86 kpc) to 0.30-Rd (= 0.68 kpc), which 
reflects that dynamical friction of SCs can cause orbital de- 
cay of massive SCs within the disk. The total number of 
SCs within 0.17?d in the disk can gradually increase so that a 
multiple SC system with the total mass (m sc ) of 4.3 x 1O 6 M 
can form in the nuclear region of the disk. The x—, y—, and 
2— components of velocity dispersions of the multiple SCs 
are 17.7 km s~ , 23.6 km s" 1 , and 3.2 km s _1 . The veloc- 
ity dispersions comparable to internal velocity dispersions 
of their nuclear massive SCs suggest that the nuclear SCs 
can soon merge to form a single nucleus. 

Figs. 2 and 3 show orbital evolution projected onto the 
x-y and x-z planes for SCs with two different masses. In both 
case, the distances from the center of the stellar disk (R) 
become significantly smaller owing to orbital decay caused 
by dynamical friction against disk field stars. The extent to 
which a SC can decay its orbit owing to dynamical friction 
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depends on its initial mass and location with respect to the 
center of its host galaxy: the SCs with masses of 3.3 x 1O 5 M0 
and 1.3 x 10 6 M Q can decrease R by 31% and 63% of their 
initial R, respectively. The orbit of the more massive SC, 
which is initially in the inner region of the disk, can be in- 
fluenced also by the inner spirals and bar developed during 
the disk evolution so that the orbital evolution can become 
more chaotic. These orbital decay seen in these two SCs can 
be more clearly seen in more massive SCs owing to more 
effective dynamical friction for them (see the Appendix A 
for details of other SCs' orbits in the standard model). 

Fig. 4 shows that about 60 % of the SCs have final 
R (Rf) significantly smaller than initial R (Ri) owing to 
dynamical friction of the SC against disk field stars. Spiral 
arms and bar in the disk can drive the angular momentum 
redistribution of the disk so that some SCs have Ri < Ri: 
these are transferred to the outer regions of the disk by 
global galactic dynamics. Fig. 4 also shows that more massive 
SCs with m sc ^ 2 x 10 5 Mq are more likely to have smaller 
Rf/Ri (i.e., more significant orbital decay due to dynamical 
friction). The mean Rf/Ri for the 20 SCs at T = 0.64 Gyr 
in this model is ~ 0.85. 

Fig. 5 shows that all SCs can significantly decrease their 
Rmin (corresponding to apo-center distances from the galaxy 
center) owing to (i) dynamical friction and (ii) radial mo- 
tions caused by non-axisymmetric structures (bar and spiral 
arms). Fig. 5 also shows that SCs with smaller Ri can expe- 
rience more significant orbital decay (i.e., smaller R m in/ Ri) 
due to dynamical friction, because dynamical friction is more 
effective in higher densities of background field stars. The 
mean R ml n/Ri for the 20 SCs at T = 0.64 Gyr in this model 
is ~ 0.34. 

Fig. 6 shows the time evolution of the total mass of SCs 
(in units of M d ) within 0.1 x Rd, which is referred to as M nuc 
for convenience. This M nuc does not literally correspond to 
the mass of a single stellar nucleus and thus goes up and 
down rapidly according to time evolution of their orbits. 
Clearly, around T = 0.4 Gyr, M nuc /M d can exceed the value 
of 0.003 observed for nucleated dwarf galaxies in the Virgo 
cluster (Cote et al. 2006). This relatively rapid formation 
of a multiple SC system is due to very efficient dynamical 
friction of massive SCs with the mean mass of 1.1 x 10 Mq in 
the disk. Therefore, the time evolution of M nuc /Md depends 
strongly on the adopted initial mass function of SCs (this 
point is discussed later). 

3.2 Parameter dependences 

Dynamical friction processes of SCs in disk galaxies de- 
pend strongly on model parameters such as Md, Mdm/Md 
(— l//d — 1), M^/Md (=/b), central stellar densities of disks, 
and M v , high- The derived M d -dependence is the most impor- 
tant result in the present study so that it is briefly summa- 
rized in the Table 2. We illustrate the derived dependences 
on model parameters as follows: 

(i) Orbital decay of SCs due to dynamical friction is 
more important in less massive host galaxies, mainly be- 
cause relative velocities (v le \) between SCs and disk field 
stars (which are comparable to stellar velocity dispersions of 
the disks) are smaller in the galaxies: the frictional drag force 
can be proportional to v rci (see 7.18 in Binney & Tremaine 
1987). Fig. 7 shows that Rf/Ri is systematically smaller in 



disk galaxies with smaller Md for a fixed fd- Fig. 7 also shows 
that for disks with Md ^ 5x 1O 1O M0, no/few SCs can show 
significant orbital decay (Rf/Ri < 0.4), which means that 
dynamical friction of SCs is not so important in their or- 
bital evolution. 

(ii) Dynamical friction of SCs by disk field stars is less 
effective in galaxies with smaller / d , mainly because stel- 
lar densities of disk field stars are lower so that frictional 
drag force can be weaker. Fig. 8 shows that for disks with 
Md/Md = 49, only a few SCs can show significant orbital 
decay (Rf/Ri < 0.4), which implies that dynamical friction 
of SCs can not be so important in orbital evolution of SCs 
in the early phase of disk evolution of galaxies when the 
mass fractions of disk stars can be very small. Fig. 8 also 
shows that a smaller fraction of SCs have Rf/Ri > 1, which 
reflects that radial redistribution of SCs due to the presence 
of non-axisymmetric structures (spirals and bars) does not 
happen in galaxies with smaller degrees of self-gravitating. 

(iii) Galactic bulges appear to be less important in or- 
bital decay of SCs within disks in comparison with other 
galactic properties such as Md and fd- Fig. 8 shows no 
remarkable differences between models with different /b 
(=0.17 and 1.0), though /b can be an important parameter 
which determines whether stellar bars can form and thus 
affect orbital evolution of SCs. 

(iv) The time scale of dynamical friction of SCs by disk 
field stars is longer in LSBs than in HSBs, because stellar 
densities of background field stars can be key factors for 
the effectiveness of dynamical friction. As a result of this, 
formation of stellar nuclei via SC merging can proceed more 
slowly in LSBs. For example, M nuc can become ~ 1O 6 M 
only after 2.4Gyr orbital evolution of SCs in the LSB model 
in which Rd is by a factor of 2.5 larger than that in the 
standard model (and other model parameters are the same 
as those in the standard model). 

(v) Final properties of SCSs (e.g., Rh. sc ) and the for- 
mation processes of nuclear multiple SCs (thus nucleus for- 
mation processes) in galaxies depend on initial MFs (e.g., 
My, high)- For example, Rh,sc can evolve more rapidly and 
have smaller final values in models with smaller M v ,hi g h (i.e., 
larger m sc for the most massive SCs) in the present study 
(e.g., i?h,sc = 0.2i? d for models with M v ,hi g h = — 11 mag). 
Also final M nuc is larger in models with smaller Mv.mm (e.g., 
M h ,nuc/M d = 0.02 for M v ,hi g h = -11 mag). This result im- 
plies that formation processes of stellar galactic nuclei (e.g., 
whether dwarfs are nucleated or non-nucleated) depend on 
original MFs of SCs in galaxies. 



4 DISCUSSION 

4.1 Evolution of SC luminosity functions 

The present study has first shown that more massive SCs 
can spiral-in to the nuclear regions of galactic disks owing to 
dynamical friction within less than ~ 1 Gyr in low-mass disk 
galaxies. It is highly likely that these multiple SCs in nuclear 
regions of disk galaxies can merge with one another to form 
stellar galactic nuclei (or single massive nuclear SCs). If this 
is the case, then the merged SCs (i.e., stellar galactic nuclei 
or nuclear SCs) would not be identified as ordinary SCs 
in observations. Therefore, it is meaningful for the present 
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Figure 9. The initial (blue) and final (red) distributions of m BC 
in the standard model. The details of the method to derive these 
SCMFs are given in the main text. 



0.006 



0.004 



0.002 



J I I I L_ 

9 

log 10 M„ (M e ) 



■ 

10 



Figure 10. The dependence of M nuc /M^ on M<j for five mod- 
els with different M d . The dotted horizontal line indicates the 
observed value (= 0.003). 



study to investigate SCMFs located within R 0.1/Zd for 
comparing the simulated SCMFs with the observed ones. 
In the followings, we assume that MFs are for SCs located 
within R ^ 0.1i?d (i-e., excluding the nuclear SCs in the 
simulations) . 

Fig. 9 shows that no SCs with log lo (m ac /M0) > 5.9 
can be seen in the final MF for the standard model, be- 
cause more massive SCs spiral-in to nuclear regions of disk 
galaxies more rapidly owing to dynamical friction of SCs 
against disk field stars. This result implies that MFs of SCs 
for 5.0 < log lo (m sc /M0) (i.e., higher-mass end) can steepen 
significantly in comparison of their original MFs (at the 
epoch of SC formation) owing to dynamical friction of SCs. 
Furthermore, given that dynamical friction of SCs against 
disk field stars is more effective in less massive disk galaxies, 
the MFs can more significantly steepen in less massive disk 
galaxies. Thus, it is highly likely that the SCMFs can be 
steeper in less massive disk galaxies. 

Based on observational studies of GCLFs for the early- 
type galaxies in the ACS Virgo Cluster Survey, Jordan et al. 
(2006) showed that the GCMFs for relatively high masses 
(m sc ^ 3 x 1O 5 M ) are steeper in less luminous galaxies 
in the Virgo cluster. They suggested that since GCMFs are 
not strongly affected by dynamical evolution of GCs (e.g., 
dynamical friction) over a Hubble time, the steeper GCMFs 
arise from variations in initial conditions. Our results how- 
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Figure 11. The same as Fig. 4 but for halo SCs in the standard 
model. Since 9 out of 20 have Rf/R[ > 1.2, they can not be seen 
in the plot. Since i?j of a halo SC is not set to be the pericenter 
distance of its orbit initially, R[ can be significantly larger than 
Ri during its orbital evolution. The mean value of Rf/R\ for all 
SCs is 1.5 and that for SCs with Rf/Ri 1 is 0.64 in this model. 
These suggest that dynamical friction of halo SCs is less effective 
than that of disk SCs in disk galaxies. 



ever imply that dynamical friction of SCs against disk field 
stars is so effective as to change the shapes of GCMFs within 
a relatively short time scale. Our results also imply that 
if these less luminous early-type galaxies were formed from 
disk galaxies either by galaxy interactions or by major merg- 
ers, their GCMFs can be steeper owing to changes of GCMFs 
in their progenitor disks caused by dynamical friction. 



4.2 Formation of stellar galactic nuclei 

Since Tremaine et al. (1975) first investigated dynamical 
evolution of GCs interacting with background stars through 
dynamical friction in M31 to discuss nucleus formation in 
M31, many authors discussed formation of stellar galactic 
nuclei via mergers of SCs (or GCs) in galaxies with different 
physical properties (e.g., Lotz et al. 2001; Capuzzo-Dolcetta 
& Vicari 2005). Recent fully self-consistent N-body simula- 
tions have shown that merging of SCs in the inner regions 
of galaxies lead to the formation of stellar galactic nuclei 
with their projected radial density profiles similar to the 
observed ones (e.g., Capuzzo-Dolcetta & Miocchi 2008). Al- 
though these works did not discuss potential importance of 
dissipative gas dynamics and star formation in the forma- 
tion of stellar galactic nuclei (e.g., Bekki et al. 2006; Bekki 
2007), they significantly advanced our understanding dy- 
namical evolution of multiple SC systems in galaxies. 

One of remaining questions in the SC merger scenario 
of nucleus formation is whether the scenario can explain the 
observed relation between stellar nucleus masses and their 
host galaxy ones (M nuc ~ 0.003M g , where M g is the galaxy 
mass; Cote et al. 2006). Since the observational study by 
Cote et al. (2006) derived total mass of the luminous com- 
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ponents (i.e., not those for entire galaxies including the ex- 
tended dark matter halos) of galaxies by using reasonable 
M/L, it is reasonable to use Md rather than M ga i for the 
purpose of discussing the observed M nuc = 0.003M g rela- 
tion. Fig. 10 shows M nU c/Md as a function of Md for the 
five models in which only Md are different between them. 
The simulated M nuc /Md is close to the observed value of 
0.003 for less luminous disk galaxies with Md < 1O 9 M . 
SCs can lose a significant fraction of their original masses 
during their orbital evolution in galaxies (e.g., Fujii et al. 
2008), which is not considered at all in estimating M nuc /M d 
in the present study. Thus the simulated M nuc /Md for less 
luminous disks can be regarded as fairly consistent with ob- 
servations. 

Fig. 10 also shows that formation of stellar galactic nu- 
clei with Mnuc/Md ~ 0.003 by SC merging is not possi- 
ble in more luminous disks owing to ineffective dynamical 
friction. Therefore, the result implies that if the observed 
M nuc = 0.003M g relation can be hold for all galaxies, the 
SC merger scenario may well have a serious problem in ex- 
plaining the observation. However, galaxies are more likely 
to have massive black holes (MBHs) rather than stellar nu- 
clei in more luminous galaxies with M g ^ 1O 1O M0 (e.g., Fer- 
rarese et al. 2006), though it is not so clear whether these 
galaxies with MBHs also have stellar nuclei (i.e., MBHs and 
stellar nuclei coexist in galaxies). Thus it would be currently 
reasonable to claim that the SC merger scenario is a promis- 
ing one for the formation of stellar galactic nuclei in galaxies. 

4.3 Compassion with dynamical friction by 
galactic halos 

We have first shown that dynamical friction of SCs against 
disk field stars can be so effective in low-mass disk galax- 
ies so that orbital decay of SCs can be as significant as 
to change physical properties (e.g., half- number radii and 
SCMFs) of their SCSs. Previous theoretical and numerical 
works focused exclusively on dynamical friction of halo SCs 
(e.g., old GCs) against dark matter halos and showed that 
the time scale of orbital decay of SCs due to dynamical fric- 
tion against dark matter halos is relatively long (a few to 
several Gyr) even for low-mass dwarfs like the Fornax (e.g., 
Goerdt et al. 2006 for the cuspy NFW profile of its dark 
matter halo). Since they did not investigate the influences 
of stellar disks on orbital evolution of halo SCs through dy- 
namical friction, it is worthwhile for the present study to 
investigate how halo SCs evolve in disk galaxies if they can 
experience dynamical friction against both dark matter ha- 
los and against disk field stars. 

Fig. 11 shows that orbital decay of halo SCs due to 
dynamical friction against halos and disks is much less sig- 
nificant owing to (i) lower mass densities of the halos and (ii) 
larger velocity differences in the SCs and dark matter halos 
for the standard model. The mean R{/Ri for the halo SCs is 
1.53, which is by a factor of ~ 2 larger than that for the disk 
SCs. Only some more massive SCs with m sc ) 2 x 1O 5 M0 
can show significant orbital decay (i.e., decreasing apocen- 
ter distances of their orbits) for halo SCs. This less effective 
dynamical friction in halo SCs is a common feature in the 
present models, which suggests that the presence of disks 
in galaxies does not significantly change dynamical friction 
processes of halo SCs. It should be however stressed that the 



halo SCs initially in the inner region of the disk galaxy in 
the standard model shows a slightly flattened final spatial 
distribution: this can be due to relatively stong dynamical 
friction between halo SCs and disk field stars in the inner 
regions of galaxies. 



5 CONCLUSIONS 

We have investigated orbital evolution of SCs being influ- 
enced by dynamical friction of disk field stars based on 
galaxy-scale numerical simulations on dynamical evolution 
of disk galaxies. We have investigated models with variously 
different model parameters in order to understand depen- 
dences of the effectiveness of dynamical friction of SCs on 
physical properties of disk galaxies. We summarize our prin- 
ciple results as follows. 

(1) Dynamical friction of SCs against disk field stars is 
much more effective in orbital decay of SCs in comparison 
with that against galactic halos in disk galaxies. For exam- 
ple, the half-number radius of the SCS (-Rh.se) in a galaxy 
with M d = 1O 9 M can decrease by ~ 30% owing to or- 
bital decay of SCs caused by dynamical friction within well 
less than 10 9 yr for a reasonable MF of SCs. However or- 
bital decay of SCs due to dynamical friction can be more 
clearly seen only for more massive SCs (e.g., those with 
m sc 2 x 10 5 M Q ). 

(2) Dynamical friction processes of SCs against disk 
field stars depend on physical properties of their host galax- 
ies. For example, the dynamical friction is much more effec- 
tive in disks with lower Md owing to smaller stellar veloc- 
ity dispersions. As a result of this, Rh, S c can more rapidly 
evolve with time in disks with lower Md . Evolution of SCSs 
in more massive disks (with Md ~ 10 10 M Q ) can not be 
clearly seen in the present study, which implies dynamical 
friction is important only for evolution of SCSs (e.g., iih,sc) 
in less massive disk galaxies. 

(3) Dynamical friction of SCs by disk field stars is more 
effective in disks with higher degrees of self-gravitating of 
the disks (e.g., higher /a), which implies that higher stellar 
densities and non-axisymmetric structures such as bars and 
spiral arms can promote dynamical friction. Orbital decay 
of SCs due to dynamical friction in disks does not depend 
so strongly on bulge masses of the disks for < /b < 1 in 
the present study. Evolution of SCSs due to orbital decay 
of SCs caused by dynamical friction is more important in 
HSBs than in LSBs. 

(4) More massive SCs can spiral-in to the nuclear re- 
gions (R < 0.17?d) of disk galaxies owing to dynamical fric- 
tion so that they can form multiple SC systems there in 
less massive disk galaxies with Md ^ 1O 9 M0. These sys- 
tems show velocity dispersions comparable to or less than 
internal stellar velocity dispersions of SCs so that they are 
likely to merge quickly with one another to form single mas- 
sive SCs there. Thus formation of stellar galactic nuclei via 
merging of more massive SCs is possible in less massive disk 
galaxies. 

(5) The observed M nuc — M g relation (i.e., M nuc ~ 
0.003M g ) for less luminous galaxies can be closely associ- 
ated with formation of stellar galactic nuclei by SC merg- 
ing. The LF/MFs for more massive SCs in SCSs of disk 
galaxies can steepen owing to the formation of stellar galac- 
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tic nuclei through merging of more massive SCs inwardly 
transferred to the nuclear region of galaxies. This steeping 
of SCLF/MFs due to dynamical friction can be more signif- 
icant in less massive galaxies. 
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Figure Al. Time evolution of R/Rd of five SCs selected from 
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m sc ^ 2 X 10 5 Mq (upper) in the standard model. Owing to dy- 
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SC-SC mutual interaction, R/Rd of these SCs can evolve in a 
complicated way. 
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APPENDIX A: ORBITAL EVOLUTION OF 
INDIVIDUAL SCS 

Fig. Al shows time evolution of distances (R in units of 
Rd) of ten SCs from the center of their host disk galaxy in 
the standard model. In order to show more clearly orbital 
evolution of each individual SC, only 5 SCs are selected and 
shown for less massive SCs (m sc < 2 x 10 5 M Q ) and for more 
massive ones (m sc > 2 x 1O 5 M ). The orbital eccentricity of 
a SC in the present study is described as follows: 



Corb — 



Rn 



Rn 



(Al) 



Rmux + R mm 

where i? m i n and i? max are the minimum and maximum R 
for the SC. It is found that i? m in of the selected ten SCs 
range from 0.006 to 0.51 and the mean R m in is 0.17. It is 
also found that e orb of the SCs range from 0.19 to 0.97 and 
the mean e or b is 0.67. 



Owing to dynamical interaction between the non- 
axisymmetric structures (e.g., bar and spiral arms) devel- 
oped in the later phase of the disk's dynamical evolution, 
orbital evolution of these SCs can be quite complicated. One 
example of this is that not all of the SCs show significant 
declines in 7? owing to dynamical friction within ~ 0.6 Gyr. 
Some of SCs show orbital decay before the formation of a 
stellar bar, their R can again become larger after the bar 
formation owing to redistribution of angular momentum by 
the bar. It appears that SC-SC interaction does not play a 
vital role in the orbital evolution of SCs owing to the small 
number (= 20) of SCs in this model. 
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